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Abstract. We investigate the heavy quark photoproduction based on the k±- 
factorization approach, focusing on the results from the saturation model. The 
deviations in the results using the unintegrated gluon distribution considering the 
saturation model and the derivative of the collinear gluon distribution are analysed. 
Total cross sections and px distributions are analysed in detail, setting the deviations 
between the color dipole approximation and the complete semihard approach. 



1. Introduction 

The investigation of heavy-quark production in high energies provides a better 
understanding of the hadron internal structure. In particular, the heavy quark masses 
are large enough to be taken as a hard scale, making the strong coupling constant 
small and allowing a perturbative QCD treatment of the process. Heavy quarks are 
produced in the clean 7p — > QQ reaction at HERA, where the incident photon (real 
or virtual) probes the proton target at high center of mass energy W . There is a large 
amount of data on heavy flavor production at HERA, which have been a plethora on 
studies analyzing the interface between hard and soft regimes jQ |2|. However, data on 
open heavy quark production are limited by small statistics and the intermediate energy 
interval between the fixed target experiments and the high energy at HERA have not 
been covered yet. In hadroproduction at the pp Tevatron collider, the situation is 
slightly better, having precise measurements of the transverse momentum distribution 
of the produced heavy quarks ||. [| . 

The underlying mechanism for heavy quark production at HERA is photon gluon 
fusion: a photon coupled to the scattered electron interacts with a gluon from the 
proton by producing a quark antiquark pair, e.g. cc (charm) or bb (bottom). In the 
collinear QCD approach, based on the well known collinear factorization theorem |5|], 
the process is described through the convolution of on-shell matrix elements, encoding 
the partonic subprocesses, with the parton distribution functions. At high energies, the 
gluon is the parton which drives the dynamics, evoluting on the virtualities obeying 
the DGLAP evolution equations M. The transverse momenta of the incident particles 
are taken as zero and in the computation of the cross sections one averages over two 
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transverse polarizations of the initial gluon. A computation of the process requires the 
knowledge of the gluon momentum distribution in the proton and the calculation of the 
7# subprocess. The gluon density as a function of the longitudinal momentum fraction x 
and virtuality scale Q 2 is known to an accuracy of a few percent from a global analysis of 
scaling violations of the proton structure function F 2 measured at HERA |7J, |], |9| . The 
QCD matrix elements of the partonic process have been calculated in next-to-leading 
order (NLO) accuracy ||10|| . Still substantial theoretical uncertainties come from the 
heavy quark mass and the from renormalization and factorization scales, hr and \xf- 

The collinear factorization approach has produced a successful description of single- 
particle distributions and total cross sections for heavy-quark production. On the other 
hand, despite of many phenomenological successes, results within this approach are in 
contradiction with data on azimuthal correlations and on distributions in transverse 
momentum of the produced heavy-quark pair |l]J (for quarkonium production see e.g. 



Refs. 0, 13] )• This problem is in general cured by introducing an intrinsic transverse 
momenta distribution (intrinsic k±) for the incident partons, parametrized using a 
gaussian profile. However, the mean value of intrinsic k± needed to describe azimuthal 
correlations and pt spectrum can be as high as (k±) ~ 1 GeV or even 2 — 3 GeV, which 
is not suitable from non-perturbative arguments. Moreover, recent calculations using 
known NLO results |13| of bottom hadropro duct ion in the collinear approach undershoot 
Tevatron data |TJ]] by a factor 2 or 3, suggesting that an important contribution to the 
computed observables is missing. 

At high energies another factorization theorem emerges, the fc^-factorization or 
semihard approach [T5| , [TE| , [17]]. The relevant diagrams are considered with the 
virtualities and polarizations of the initial partons, taking into account the transverse 
momenta qi± and q 2 i_ of the incident partons. The processes are described through 
the convolution of off-shell matrix elements with the unintegrated parton distribution, 
T{x 1 k±). The latter can recover the usual parton distributions in double logarithmic 
limit (DLL) by its integration over the transverse momentum k±. In the asymptotic 
energy limit, the unintegrated gluon distribution should obey the BFKL evolution 



equation [ 18[ . At the present, there is a lack of an accurate determination of this 
quantity at the same level as the usual parton densities. The matrix elements computed 
for the relevant subprocesses within the fc^-factorization approach are more involved 
than those needed in the collinear approach already at LO level. On the other hand, a 
significant part of the NLO and some of the NNLO corrections to the LO contributions 
on the collinear approach, related to the contribution of non-zero transverse momenta 
of the incident partons, are already included in the LO contribution within the k±- 
factorization formalism [0]. Moreover, part of the virtual corrections can be resummed 
in the unintegrated gluon function |TI|] . Recently, the computation of NLO corrections to 
the subprocesses have been completed and full calculations at this level should appear 
in the near future [Tj|. An important feature of the approach is the equivalence at 



leading logarithmic approximation with the color dipole picture P0| , which have been 
used in a profuse phenomenology at HERA Besides, a very important issue is 
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the consistency of the fcj_-factorization approach including nonleading-log effects with 
the collinear factorization beyond leading order p^[ : the coeficient functions and the 
splitting functions providing q(x, Q 2 ) and G(x, Q 2 ) are supplemented with the all- 
order resummation of the a s ln(l/x) contributions at high energies, in contrast with 
a calculation in fixed order perturbation theory. 

Two additional ingredients should be taken into account when considering the 
semihard approach: the infrared sector and saturation effects. The unintegrated gluon 
function should evolve in tranverse momentum through the BFKL evolution at high 
energies, leading to the diffusion on k^ of the initial gluons in the evolution process 
|23|. In this diffusion scenario the transverse momenta values are spread out into the 
infrared (and ultraviolet) region, where the perturbative description is not completely 
reliable. Therefore, the evolution should take into account properly the correct behavior 
in that region. The recently calculated non-linear corrections to the BFKL approach 
[ 24| introduce a natural treatment for these difficulties, where the saturation scale Q s 
provides the suitable cut-off controlling the infrared problems. As the longitudinal 
momentum fraction x ~ Q 2 /W 2 decreases, unitarity corrections become important and 
control the steep growth of the gluon distribution. In this domain the gluon distribution 
can saturate completely or acquire a mild logarithmic increasing. The most appealing 
approach taking into account both the notions of infrared behavior (confinement) and 
parton saturation phenomenon is the saturation model which is an eikonal-type 
model based on the color dipole picture of high energy interactions. In this approach 
the physical picture is quite simplified and the expression for the saturation scale is 
promptly calculated. The adjustable parameters of the model are obtained from a fit to 
small-x HERA data on the inclusive structure function and the photoproduction total 
cross section, being suitable for further applications to more exclusive quantities as open 
heavy-quark production. 

This paper is organized as follows. In the next section we present the main formulae 
for the heavy-quark photoproduction within the fc^-factorization approach, defining 
the relevant kinematical variables. Furthermore, we investigate the unintegrated gluon 
distribution from the saturation model in comparison with the results from the derivative 
of the integrated gluon distribution, xQ(x, k±). In Sec. (||), we present our results on 
the total cross section for charm and bottom photoproduction, as well as estimates 
for the p^-distribution of the produced heavy quark pair. We investigate in detail the 
deviations in the results when one confronts the dipole approximation and a conservative 
fe^-factorization procedure. Namely, the choice of scale for the coupling constant and the 
suitable longitudinal momentum fraction entering in the unintegrated gluon function. 
The predictions from the derivative of the collinear gluon function are also studied. In 
the last section, one presents the conclusions and final considerations. 
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2. Heavy-quark photoproduction in the fc^-factorization 

In this section we investigate the quasi-real photon scattering off a proton in the 
semihard approach. In light-quark photoproduction there is a poor knowledge 
concerning the effective light quark mass, which is associated with non-perturbative 
aspects of the process. This problem is naturally solved in heavy-quark photoproduction, 
because of the heavy quark mass. The semihard approach is valid in the domain 
where the following double inequality holds: s ^> fi 2 ~ s ^> Aq CD , i.e. the 
typical parton interaction scale /i 2 is much higher than the QCD cutoff parameter 
Aq CD , and much lower than the center of mass energy, ^/s. The fc^-factorization 
approach can resum in the leading logarithmic approximation all the contributions 
proportional to [a s \n(fi 2 /AQ CD )] n , [a s \n(n 2 / Aq CD ) \ia(l/x)} n and [a s ln(l/x)] n , where 
the first one corresponds to the collinear DGLAP resummation 0, the second one the 
double logarithmic contribution and the last one the BFKL resummation [Q. Such 
resummation leads to the unintegrated gluon distribution, J-{x, k±), which can also 
depend on the scale /i 2 . It gives the probability to find a parton carrying a longitudinal 
momentum fraction x and transverse momentum k±. Requiring that in the DLL 
limit, the transverse photoabsortion total cross section written in the fcj_-factorization 
approach to be consistent with the same limit from the DGLAP approach, the collinear 
gluon distribution is given by the unintegrated one in the following way, 

xG(x^ 2 ) = J ^F( x ,k 2 ± ), (1) 

where /i 2 is the scale of the process, which can be Q 2 , the heavy-quark mass or the p± 
of the produced particles, for instance. In particular, in obtaining Eq. ([I]), a strong 
ordering condition k 2 ± <C p\ <C /i 2 is considered, where p ± is the transverse momentum 
of the quark-antiquark loop in the photon vertex. 

In order to compute the cross section of a physical process |27|] , the unintegrated 
gluon function should be convoluted with the off-shell matrix elements for the relevant 
partonic subprocesses. In these matrix elements the polarization tensor of the virtual 
gluon is given by Lffi = Eg e* u = k^_k u ± /\kj_\ 2 . In the following we will calculate the 
total and differential cross section of heavy-quark photoproduction (charm and bottom) 
taking into account the diagrams shown in Fig. ([j]). For convenience, one can define 
the Sudakov variables for the process ep — > QQX at high energies, 

pi = a 1 P 1 + j3xP 2 + p 1± , q = X\P\ + q ± , (2) 

p 2 = a 2 Pi + P2P2 + P 2 ± , k = x 2 P 2 + k ± , (3) 

where, as it is shown in Fig. ([!]), p\ and p 2 are the four-momenta of the produced 
heavy-quarks, q and k are the photon and the gluon four-momenta, respectively. The 
corresponding transverse momenta are p 2 j_, q± and k±. The electron and proton 
momenta are denoted by Pi and P 2 . In the center of mass frame of the process, one has 
P l = (VS/2, 0, 0, P 2 (v^/2, 0, 0, -vty2), P 2 = Pf = and (P x • P 2 ) = s/2, with 

yfs being the center of mass energy. From simple inspection of Eqs. @ and @ and 
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Figure 1. The leading order QCD diagrams for heavy-quark production via photon- 
gluon fusion. The momenta of the particles are shown and the blob represents the 
gluon emission chain encoded in the unintegrated gluon distribution, J-{x2, kj_). 



using the energy-momentum conservation law, one obtains the relations, 



q- = q- ± , K- = k{, (4) 
q± + k ± = p 1± + p 2± . (5) 

The Sudakov variables can be written in terms of the transverse masses m\ 2 ± = 
mn +Pi2j_, where niq is the heavy-quark mass, and the heavy-quark rapidities y* 2 , in 



Q 



the following way 
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where in the photoproduction case q = Pi, the variable x\ simplifies to x\ = 1 and 
one can define ol\ = z and a 2 = 1 — z. The z and (1 — z) variables correspond 
to the longitudinal momentum fraction carried by the heavy-quark having transverse 
momentum p l± and p 2J _, respectively. 

Having introduced the relevant definitions and variables, the differential cross 
section for the photoproduction process is expressed as the convolution of the 
unintegrated gluon function with the off-shell matrix elements ITJI, |2( 



da{ 1P -> QQX) 



d 2 p 



i± 



dy\ d 2 k 



F{x 2 ,k 2 ) I M 1 2 (off -shell) 



ira 2 



(9) 



where the off-shell LO matrix elements are given by J28, p6 



|A^| 2 (off — shell) = a em a s (fi 2 )e 



Q 



z 2 + (l-z) 



[(P2 - k) 2 - m 2 Q ] [(pi - k) 2 - m 2 Q ] 
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+ 



2m 2 Q ( : 1 - : 



k\ \ (pi - k) 2 - m 2 Q (p 2 - k) 2 - m 2 Q 



, (10) 



where a em = 1/137 is the electromagnetic coupling constant and cq is the electric 
charge of the produced heavy-quark. The scale /i in the strong coupling constant will 
be specified latter on. In general, it is taken to be equal to the gluon virtuality, /i 2 = k 2 
in close connection with the BLM scheme [29|. In the leading ln(l/x) approximation, 



a s should take a constant value. When the transverse momenta of the incident partons 
are sufficiently smaller than those from the produced heavy-quarks, the result from the 
collinear approach is recovered. The final expression for the photoproduction total cross 
section considering the direct component of the photon can be written as 

, --) 2 ] 

where D x = p\ ± + m 2 Q and D 2 = (k± — p 1± ) 2 + ml 

In Eq. ( |Tl~D the unintegrated gluon function was allowed to depend also on the 
scale /i 2 , since some parametrizations take this scale into account in the computation 





of that quantity (see, for instance Ref. |30| for a compilation of a number of them). 



We are now ready to calculate the total and differential cross sections for the process, 
provided a suitable input for the function J r (x 2 , k 2 ± ; fi 2 ). The practical procedure in this 
paper is to consider one of the simplest parametrizations available, covering a consistent 
treatment of the infrared region and taking into account the expected saturation effects 
at high energies. These features are fullfilled in the saturation model, which it is shortly 
reviewed in the following and it will be conveniently contrasted with the results coming 
from the derivative of the collinear gluon function, 

d [xG(x, k 2 ±) 
d \nk\ 

where G(x, fi 2 ) is the integrated gluon distribution, which can be taken from the available 
parameterizations in the literature 0, || [J. 

2.1. The saturation model 



The saturation model [[25fl is based on the color dipole picture of the interaction, 
represented in the target rest frame where the transverse size r of the dipole quark- 
antiquark coming from the virtual photon Fock state fluctuations is fixed during the 
interaction. This representation can be recovered from the fc^-factorization approach 
in the leading logarithmic approximation through the Fourier transform between the 
transverse size and the transverse momentum spaces. The photoabsortion total cross 
section is written as a convolution of the virtual photon wavefunction with the effective 
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cross section for the interaction dipole-target, 



7 P 
tot 



>,Q 2 ) = Jdz Jd 2 r (\V T (z,r,Q 2 )\ 2 +\* L (z,r,Q 2 )\ 2 ) a dvp (x } r) } (13) 

where the dipole cross section interpolates between the color transparency behavior for 
small size dipole configuration and the confinement features for large dipole sizes. It 
reads 



a d 



ip 



R 2 (x) 



0-Q 



1 — exp 



AR 2 (x) 



x 

•i'O 



(14) 



(15) 



GeV 2 

where Rq{x) is the saturation radius which decreases when x —>■ 0. The parameters ao, 
xq and A are determined from a fit to small-x HERA data. An additional parameter 
is the effective light quark mass m q = 140 MeV, which is needed to produce finite 
results on the photoproduction total cross section. The saturation scale is defined as 
Q 2 = 1/-Rq(x): when rQ s /2 <C 1 the model reproduces color transparency, o~di V ~ r 2 , 
whereas when r Q s /2 > 1 the cross section tends to a constant at large r Q s (it simulates 
confinement), a dip ~ a . In the region rQ s /2 ~ 1, the model simulates the physics 
from the multiple scattering resummation of gluon exchanges in an eikonal-type way 
representing the black disk limit of the proton. In the original model, the coupling 
constant is considered fixed as a s = 0.2 and in order to consider the formal limit of 
photoproduction, the Bjorken variable has been modified in the following way, 



x 



X B j 1 + 



Am' 



Q 2 + Am 



(16) 



Q 2 J W 2 

Starting from the fc^-factorization approach for the total photoabsortion cross 
section, but disregarding a transverse momentum dependence in the argument of the 
strong coupling constant and in the variable x, an analytical result for the Fourier 
transform between the spaces can be obtained, 

2tt f d 2 k_ 



o- d ip{x, r) 



a s J-"(x, k 2 



1-e 



ik 



1-e 



ik 



(17) 



Therefore, the Eq. (\l7\j can be used to extract the unintegrated gluon function 
from the model in the r-space, once the dipole cross section has a finite limit at r — > 0, 



denoted as <J^ p \x) 



a s ^(x, k 



It can be written as H31 
d 2 r 



3 
An 
3 



exp (i k± ■ r) 



(2vr) 2 

drr Jo(k± r 



k 



k\ 



(19) 



»^ JO 

After writting down the expression for the unintegrated gluon function to be 
employed in the studies of the next section, some considerations are in order. The 
saturation model is a quite successful approach for the small-x region and it was also 
extended to simultaneously describe the diffractive DIS. However, when x — > 1, the 
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approach is no longer suitable and threshold correction factors should be introduced. 
The simplest way to implement this is to consider dimensional-cutting rules: for a 
subprocess having n spec spectators quarks which do not interact with the photon, the 
corresponding threshold factor is given by (1 — x) 2nspec ~ 1 . For example, including light 
and charm quarks, the number of spectators is n spec = 4 and in our analysis we include 
the multiplicative correction factor (1 — x) 7 in the unintegrated gluon function from 
the saturation model. In next section, such a correction provides a correct description 
of the fixed target energy region, whereas the results from the original model remain 
unalterated at intermediate and high energies. 

To our purposes in this investigation, we use the following parameters corresponding 
to the parameterization which includes the charm quark: <t = 29.12 mb, A = 0.277 and 
Xq — 4.1 • 10~ 3 From Eqs. ([i~8|) and (0), the unintegrated gluon distribution from 
the original saturation model supplemented by the threshold factor is given by 

^kl) = ^^R 2 (x)k 4 L ex P (-R 2 (x)kl) (1-x) 7 . (20) 

~r 7T Gig 

We illustrate in Fig. (Q), shown in the plot on the left, the unintegrated gluon 
function from Eq. (|20D as a function of the transverse momentum k\ at typical values 
of the x variable, covering large and small longitudinal momentum fractions. We remark 
that the small- x region corresponds to x < 10~ 2 . Its main features are quite clear: the 
function is peaked at k\ = 2Q 2 , with a narrow distribution around this value, being 
somewhat slightly asymmetric for large x. Therefore, the peak is shifted to larger k\ 
as x decreases. The most important feature is a large contribution from the very small 
transverse momentum sector, k^ < 1 GeV 2 at large x. Moreover, the unintegrated gluon 
distribution is strongly suppressed for large k\ due to the missing parton evolution in 
the original model. This shortcoming is cured in the recent implementation of DGLAP 
evolution for the unintegrated gluon distribution BT . 



In order to investigate the importance of a fe^-enhancement coming from a QCD 
parton cascade emission, we present in the plot on the right in Fig. (Q) the result from 
the derivative of the collinear gluon distribution, Q{x, k 2 ± ), where it was chosen the LO 
GRV98 gluon parametrization for the calculation. The use of this quantity gives the 
possibility to check the consistency of introducing elements from DGLAP evolution in 
the fc^-factorization approach. That parametrization has the virtuality Qq = 0.8 GeV 2 
as initial evolution scale and bellow this value one should make an assumption for the 
behavior of Q. In our analysis we apply the following procedure, 



Q(x, k]_) = k 



2 d[xG(x,k 



2 



dk 



2 



e(Ql - k] 



+ d[ 7!nkl ±)] Qikl - Qlh (21) 

where the first term can recover xG(x, Q^) by integration over transverse momentum 
from up to Q^. However, we emphasize that such a procedure is not unique and 
other ansatze can be introduced. For instance, we can consider a lower integration limit 




Figure 2. Left plot: The unintegrated gluon distribution from the saturation model 
(GBW) as a function of k\ for typical x values. Right plot: the derivative of the 
collincar gluon distribution (LO GRV98 parametrization) as a function of k\ for typical 
x values. 



k±min — Qli emphasizing that in the theoretical curves the cut-off mainly would affect 
the overal magnitude of the cross section. As expected, in our results the transverse 
momentum spectrum is broader, in contrast with the saturation model. At the small- a; 
region the deviation is huge both in behavior and in overall normalisation. The large 
x region is correctly described since the collinear gluon function is adjusted on the 
whole kinematical range at HERA. It has been verified that by using the NLO collinear 
gluon distribution the deviations are negligible f32fl . A technical remark is the pronged 



behavior of Q(x, k ± ) above Qq, which has no physical meaning, since it has to do with 
the grid interpolation routines used to obtain the collinear gluon function, and this effect 
is smoothed out in the integrated quantities and should not affect our latter results. 

In the next section we investigate the saturation model in the computation of the 
total and differential heavy-quark photoproduction cross section, and a comparison with 
the derivative of the collinear gluon function will be used in order to study the effects 
of QCD evolution. 

3. Results and discussions 



The available data [|33| [M, 35), |36| on heavy-quark (charm and bottom) photoproduction 
range from energies W of fixed target experiments about tens of GeV, up to the 
high energy HERA data around W ~ 200 GeV. Kinematically, the low energy data 
corresponds to x ~ 10 _1 and the high energy ones to x ~ 10~ 4 . The experimental 
errors are rather large and the intermediate region between low and high energies is not 
covered with measurements. Beauty production at HERA is suppressed by two orders 
of magnitude with respect to charm, due to the larger mass and smaller electric charge 
of the b quark. 
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Figure 3. Left plot: The charm photoproduction total cross section as a function of 
center of mass energy W and the fc^-factorization approach results using the saturation 
model and different technical procedures in the calculation (see text). Right plot: The 
bottom photoproduction total cross section as a function of center of mass energy W 
and the fc^-factorization approach results using the saturation model and different 
technical procedures in the calculation (see text). 



The total and differential cross section are computed from Eq. ( O ) and the 
unintegrated gluon function from the saturation model, Eq. (|20|). The heavy quark 
mass was considered as m c = 1.5 GeV for charm and mj, = 4.8 GeV for bottom. In 
order to investigate in detail the results emerging from the fc^-factorization approach 
and in an attempt to go beyond the leading logarithmic ln(l/x) approximation, we have 
considered in our analysis the following procedures: 

(i) We keep the original features of the saturation model. Namely, a fixed strong 
coupling constant a s = 0.2 and the longitudinal momentum fraction x 2 being 
the photoproduction limit of the Bjorken xsj, given by Eq. flT6|) , entering in the 
unintegrated gluon distribution. This procedure is equivalent to using the color 
dipole picture and the saturation model, as performed in the detailed study of 
Ref. ]37]. The result is shown in the dot-dashed curves in Fig. (^), for charm 



(plot on the left) and bottom (plot on the right) total cross sections. There 
is good agreement with the low energy data, whereas the high energy data are 
underestimated. The good description of the fixed target data is ensured by the 
threshold correction factor (1 — x) 7 , since the original model would overestimate 
the data already at intermediate energies. At high energies, those threshold effects 
are completely negligible at W > 20 GeV for charm and W > 50 GeV for bottom. 



Similar conclusions are obtained in the Ref. |j37|| . One small difference between 
the procedure above and the calculations in Ref. |37|] is the prescription for the 
longitudinal momentum fraction being xq = mq/z{l — z) in the latter. If we 
consider a mean value z = 0.5 for the quark momentum fraction, the results are 
completely equivalent. The very close similarity between our predictions and those 
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of Ref. |J7| corroborates our procedure of calculation. 



We allow the argument of the strong coupling constant to run with the following 
scale, /i 2 = k 2 ± + This procedure is still very close to the saturation model, 
and it is the general procedure in calculating observables by other groups [^7], [26[ . 



The term /Xq is introduced in order to avoid divergencies coming from the coupling 
constant in the infrared region. The value /i 2 , = 1 GeV 2 has been used, motivated by 
the value of the saturation scale Q s , ensuring that the low transverse momentum 
region is dominated by that scale. The result seems to spoil the previous good 
agreement at low energies, as shown in the dashed curves in Fig. ([3D. However, the 
high energy data are described in better agreement than in the previous procedure. 
In conclusion, the introduction of a running coupling in the calculation shifts the 
overall normalisation towards higher values at high energies, by enhancing the k 2 ± - 
profile in the unintegrated cross section. 

(iii) As a final aspect, we consider a conservative procedure concerning the k±- 
factorization approach. The argument of the strong coupling constant is let to run 
with the scale fi 2 = p 2 j_+rnQ, where p ± is the transverse momentum of the produced 
heavy quark. Moreover, we have a different prescription for the x<i variable than 
in Eq. (|16|) . The correct value for the longitudinal momentum fraction entering in 
J r (x2,k 2 L ) is given by the definitions in Eqs. (fD and (|7|). Namely, the momentum 
fraction is given by 

_ mf^ m 2 ± _ m 2 Q + (p ± - k ± ) 2 m 2 Q + p\ 
%2 ~ zW* + (l-z)W*~ z + l-z ' { ' 

where one has made use of the relations in Eqs. (|6|) and (0) and the energy- 
momentum conservation law, Eq. (|5]). Although Eq. (f22|) is well defined, it is 
involved for practical computations and for simplicity we rely on the following 
approximation. We benefit from the results on the fc^-factorization approach 
applied to ep collisions, in particular for the proton structure function, performed 
in Ref. f38| . There, it was verified that for not too high virtualities Q 2 (including 
photoproduction), a suitable approximation is x 2 = 2%Bj- This value is obtained 
from a careful investigation of the contributions in the transverse momentum 
integration range for the corresponding DGLAP piece and for the semihard 
approach. The results using the procedure above is shown in the solid lines in 
Fig. (^), presenting an intermediate behavior between the first and the second 
procedures. In particular, the fixed target data are described in good agreement and 
the results for high energy data are slightly improved compared to the equivalent 
dipole result (specially for the charm results). 

A more detailed calculation considering the resolved part of the photon is 
beyond the scope of the present analysis. In Ref. [27|, the resolved component is 



considered by including the off-shell matrix elements identical to those of heavy quark 
hadroproduction, convoluted with the photon and the proton unintegrated parton 
densities. In Ref. |[7|, where the color dipole picture is applied to heavy-quark 
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Figure 4. Left plot: The k± profile (charm) for the unintegrated cross section from 
the saturation model for two different procedures concerning coupling constant and 
momentum fraction scale (see text). Right plot: The k± profile (charm) for the 
unintegrated cross section from the derivative of the collinear gluon distribution (LO 
GRV98 parametrization). 



photoproduction using the saturation model, a VDM contribution is included. In 
our analysis, we only consider the matrix elements of the direct component of the 



photon. According to the authors of Ref. [39], the 7 — > cc component of the photon 



is automatically included in the fc^-factorization approach, since there is no restriction 
on the transverse momenta along the evolution chain. Moreover, it is shown in Ref. 
39| , in particular when calculating the D* (mesons) differential cross section da/dx. 
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that part of the resolved photon contribution is effectively included by the BFKL or 
CCFM [[IIJ evolutions, namely this is included in the evolution of the unintegrated gluon 
distribution, and the off-shell matrix elements would contain only the direct component 
of the photon. When using such approaches in photoproduction, the resolved photon 
contribution is in general not included explicitly in order to avoid double counting, since 
the off-shell gluon from the BFKL evolution would take already into account a certain 
portion of this contribution. 

In order to investigate the influence of the procedures above in the predictions for 
the total cross section and to find from where (in the transverse momentum range) 
comes the main contribution, we propose to consider the k^ profile of the unintegrated 
cross section. This quantity, denoted by W(x,k 2 L ), is obtained by unfolding the k\- 
integration in Eq. (|TT| ) . In Fig. (^) we show the profile functions for charm, by using the 
saturation model (plot on the left) and the derivative of the collinear gluon distribution 
(plot on the right), for fixed target energies (momentum fraction x = 10 -2 ) and high 
energies (x = 10~ 4 ). For the saturation model, we consider two of the procedures above, 
namely the dipole approximation (fixed a s and X2 = xbj) and the semihard approach 
(a s running with /z 2 = p\ + uIq and x 2 = 2xsj)- At x = 1CT 2 the profile functions 
are peaked around k 2 ± ~ 0.3 GeV 2 , whereas at x = 10" 4 the peaks are shifted towards 
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Figure 5. Left plot: The k± profile (bottom) for the unintegrated cross section from 
the saturation model for two different procedures concerning coupling constant and 
momentum fraction scale (see text). Right plot: The k± profile (bottom) for the 
unintegrated cross section from the derivative of the collinear gluon distribution (LO 
GRV98 parametrization). 

k\ ~ 1 GeV 2 . As expected from results for the total cross section (Fig. the semihard 
approach results have a larger normalization than the dipole ones. An important feature 
emerging from these results is the dominance of the small k 2 ± region for the charm total 
cross section. Indeed, at high energies the peak is of order of the saturation scale, 
k x — Q 2 , and contributions for k\ > 10 GeV 2 are negligible. 

For the derivative of the collinear gluon distribution (usually denoted d — Gluon), 
we choose the scale /z 2 = p 2 + m 2 Q and X2 = 2xsj- We can notice in Fig. (f|b) the effect 
of the discontinuity in the derivative of the gluon function, Eq. ( |2l~D at Ql = 0.8 GeV 2 : 
the profile function peaks in this value for k 2 ± > Qq, whereas it has a flatter behaviour 
for k 2 ± below Qq. We can also notice that the k± profile for the derivative of gluon 
distribution is broader than in the saturation model. Still the main contribution comes 
from the small k± region, however, intermediate values of transverse momentum give a 
considerable contribution. This will bring the predictions closer to data, as we will see 
latter on. 

In Fig. (|5[) we present the results for the k± profile for the bottom, where we 
compare the different procedures and unintegrated gluon functions, similarly to the 
charm case. Some differences in comparison with charm are evident. For the saturation 
model calculations, the results of the dipole and semihard procedure are very similar, 
since the scale /i 2 = p\ + m 2 Q provides a high virtuality even at very small transverse 
momentum due to the large mass of the bottom. This makes the strong coupling 
constant to be very close to the value a s = 0.2 in the whole p ± range, as assumed 
in the original saturation model. The results from the derivative of the collinear gluon 
distribution keep the same features as the charm case, having again a broader k± profile 
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in comparison with the saturation model. 

The study of the k^ profile shows the dominant region in the k^ range and the 
effect of choosing distinct scales for the coupling and longitudinal momentum fractions. 
From the profiles discussed above, we expect that the broader k± spectrum for the 
derivative of the collinear gluon function will enhance the overall normalisation of the 
total cross section at high energies, improving the data description in comparison with 
the saturation model predictions. Motivated by this fact, we also perform a comparison 
at total cross section level between the two unintegrated gluon functions, namely the 
saturation approach and the derivative of the collinear gluon funtion. In both cases we 
make our default choice of scale fi 2 = p\ + uiq and %2 = 2xBj- This comparison is 
shown in Fig. (|6]) for both charm and bottom total cross sections. The saturation model 
underestimates high energy data, since the treatment of QCD evolution is not considered 
in the original model. Recent improvements, taking QCD evolution into account, should 



cure this shortcoming f31| . The derivative of the collinear gluon distribution gives a 
better description of high energy data, since it includes the referred gluon emission. As 
expected, it is in disagreement at lower energies, since the non-singlet (valence) content 
was not included in the analysis. In addition, it was also verified that the corresponding 
unintegrated distribution function takes negative values in that region. For sake of 
illustration, we also show the parton model results (collinear approach) for the LO 
process jg — > QQ, where it has been used m c — 1.3 GeV, m& = 4.75 GeV, and /i 2 = s. 
This gives a reasonable description of data given the use of lower heavy quark masses 
or alternativelly considering higher order corrections to the LO calculation. In contrast, 
the semihard approach gives a reasonable description of data already at LO level. The 
energy dependence is distinct in the calculations: the saturation model provides a mild 
energy growth, whereas in the collinear approach the growth is steeper. The collinear 
approach and the semihard result using the derivative of gluon function present a similar 
energy behavior, with sizeable deviations only at low energies near the threshold. 

As a final investigation, we compute the p ± -distribution for charm and bottom using 
the saturation model (using the 3 procedures discussed earlier) and the derivative of the 
collinear gluon distribution, at center of mass energy W = 200 GeV. The predictions are 
shown in the Fig. ([?]). A remarkable feature is the finite and well controlled behavior at 
small p ± for both gluon functions. One can also see the usual fall off at large transverse 
momentum. The finiteness at zero transverse momentum in a LO level calculation is 
one of the main advantages of the semihard approach. Our results are comparable with 
those of Refs.[27, 26], which consider other parametrizations for the unintegrated gluon 
function. Our results for the saturation model are quite similar, even using different 
prescriptions for the scales at a s and for longitudinal momentum fraction, with a slightly 
deviation at larger transverse momentum. The growth at small p ± is less steep in the 
bottom case than in the charm calculation, because of the larger bottom mass in the 
argument of a s . The derivative of the collinear gluon distribution (d — Gluon) produces 
a similar behavior on p ± , but with a somewhat higher overall normalization. 



Heavy Quark Photoproduction in k ±- Factorization Approach 



15 



10 F 




W(GeV) 

Figure 6. The results for the charm and bottom total cross sections considering the 
saturation model, the derivative of the collinear gluon distribution and the collinear 
parton model. See discussion in the text. 
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Figure 7. Left plot: The charm ^-distribution from the saturation model (using 
three distinct procedures) and the derivative of the collinear gluon distribution for 
W = 200 GeV. Right plot: The bottom p ± -distribution from the saturation model 
(using 3 distinct procedures) and the derivative of the collinear gluon distribution for 
W = 200 GeV. 

4. Conclusions 

In this paper we investigate in detail the fc^-factorization approach (semihard approach) 
applied to heavy-quark photoproduction. In this formalism the cross section is given by 
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the convolution of off-shell matrix elements with the unintegrated parton distributions, 
J r (x2,k 2 L ). The matrix elements are at present known at LO accuracy, and include 
most of the NLO and even NNLO diagrams from the collinear factorization approach. 
This fact is advantageous in heavy-quark production, since NLO calculations in the 
collinear approach undershoot data (specially Tevatron data for bottom). There are 
several parameterizations for the unintegrated gluon function, relying on the solution of 
evolution equations or based on phenomenological considerations. We have investigated 
the application of the saturation model parameterization, which provide us a safe 
treatment of the infrared region and includes the onset of the parton saturation 
phenomenon. Moreover, the adjustable parameters of this model are extracted from 
the high energy HERA data, and therefore the results are parameter free. 

In order to go beyond the leading logarithmic approximation (ln(l/x)), we let 
the strong coupling constant run and use a suitable longitudinal momentum fraction 
entering on T . The description of the total cross section data is strongly dependent on 
those procedures. An additional ingredient in the calculations is a threshold correction 
factor accounting for the low energy behavior. The best agreement is obtained by the 
prescription /i 2 = p\ + m 2 Q and x<i = 2xsj- However, in any case the saturation model 
slightly underestimate the high energy experimental results. This comes from the fact 
that QCD evolution is not present in the original model. In order to investigate the role 
played by these evolution emissions, we have considered the derivative of the collinear 
gluon distribution, which provides a closer connection with the DGLAP formalism. 
Indeed, the results for the charm and bottom total cross sections are in nice agreement 
in both low and high energies. Recent improvements of the saturation model considering 
these emissions are expected to produce similar results. 

It was verified that the study of the profile provides important information where 
in the range of transverse momentum the main contribution for the processes is coming 
from. For the saturation model the most important piece is peaked at the saturation 
scale k 2 ± ~ Q 2 S . The k± spectrum is broader for the derivative of the gluon distribution, 
Q{x2i k 2 ± ). The p_|_-distribution of the produced heavy-quarks is also computed, showing 
the effects coming from different prescriptions for the scales considered. The results for 
the saturation model are quite similar, even using different prescriptions for the scales 
at a s and for the longitudinal momentum fraction, with a slight deviation at large 
transverse momentum. The growth at small p L is less steep in the botton case than for 
the charm. The derivative of the collinear gluon distribution produces a similar behavior 
on p L . 

The study of heavy quark photoproduction in the framework of the semihard 
approach improves the understanding of QCD dynamics both in the infrared and in 
the perturbative regions. It also sheds light on the proton structure, specially the gluon 
distribution. Therefore, it contributes in the understanding of the amazing interplay of 
soft and hard QCD phenomena. 
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